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Effects of externally imposed periodic changes in the environment on population dynamics are 
studied with the help of a simple model. The environmental changes are represented by the temporal 
and spatial dependence of the competition terms in a standard equation of evolution. Possible 
applications of the analysis are on the one hand to bacteria in Petri dishes and on the other to 
rodents in the context of the spread of the Hantavirus epidemic. The analysis shows that spatio- 
temporal structures emerge, with interesting features which depend on the interplay of separately 
controllable aspects of the externally imposed environmental changes. 
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INTRODUCTION 

The mathematical study of population dynamics has 
been a subject of great interest in recent years, with ap- 
plication widely spread among different fields. An exam- 
ple is the description of emerging patterns in bacterial 
colonies ll, J,, ik, ..S,, ^ JJ • Other studies describe the 
behavior of the populations of superior organisms such 
as insects and rodents j3. ifll linL . 

If, in basic models for the description of the evolution 
of a given population, we focus attention on reproduc- 
tion, competition for resources, and diffusion, the Fisher 
equation JJj, J^J appears to be a useful mathematical 
tool: 

du(x,t) ^d'^u{x,t) , . ,9, , , , 
qI ' +au{x,t)-hv?{x,t). (1) 

Diffusion with coefficient D is considered here as well as 
the growth of the population at rate a and a competition 
process weighted by 6, also containing environmental fea- 
tures. In the present paper, we consider effects of spatio- 
temporal dynamics in the nonlinear term as a represen- 
tation of externally imposed or seasonal environmental 
variations. We thus take 6 in Eq. ^ to be space and 
time dependent, b{x,t). Environmental variations could 
alternatively be considered as affecting a (such that a is 
dependent on time and space rather than b) as has been 
done in some earher studies 0, 0|- However, we restrict 
our attention to a constant a and varying b{x, t) because 
this allows us to define a Fisher velocity and size in terms 
of a constant a. Following we consider a bounded re- 
gion which will be referred to as a bubble of favorable 
conditions suitable for the survival of a given species. 
Outside of this region the environmental conditions are 
so harsh that there is no possibility of survival. 

Associated with the Fisher equation there is a nat- 
ural velocity and a natural length, 2\J Da and TTyjD/a, 
respectively. The former is known as the Fisher velocity 
and is the velocity at which fronts tend to travel after 



adjustments from most initial conditions Thus, an 
arbitrary shape of the initial population with compact 
support will have fronts that move, after initial tran- 
sients, with the Fisher velocity. It is an important quan- 
tity in experiments on bacteria with moving masks [3| as 
it represents the minimum mask velocity at which bac- 
teria tend to extinction, rather than being able to follow 
the bubble through the combined effect of growth and 
diffusion. The Fisher length tt-^/D/o arises 0, 0| in 
the bacterial context as the minimum length of the mask 
below which the steady state population of the bacteria 
vanishes as the bacteria diffuse out of the masked area 
quicker than they can grow to any finite saturation value. 
The Fisher length is thus the 'diffusion length' within the 
growth time (reciprocal of the growth rate a) . 

Consideration of these two quantities suggests that it is 
natural to envisage two observations involving externally 
imposed periodic variation of the environment. The first 
is one in which experiments of the kind reported in ref. 
Q are carried out with the velocity varying periodically, 
in other words with the bubble repeatedly returning to 
its original position. This is the oscillating bubble case. 
The second experiment is one in which the bubble extent 
is made to vary from a large size to one below the Fisher 
length (i.e. the extinction size). The objective of the 
latter (breathing bubble) experiment is to analyze ex- 
tinction tendencies under periodic variations of the bub- 
ble size. We analyze these two hypothetical experiments 
in this paper through numerical calculations based on a 
Crank-Nicholson scheme. 



TRAVELING BUBBLE 

Solutions of the Fisher equation for bounded domains, 
corresponding to static bubbles, have been already dis- 
cussed in [7] and references therein. A natural change, 
to be considered as a first approximation to the problem 
addressed here, is to allow a translation in the bubble. 
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FIG. 1: The bold curve depicts Umax clS the bubble moves 
adiabatically; acceleration is affectively zero. The remain- 
ing curve displays Umax for a bubble accelerating linearly at 
.la^Wc. 



We have performed a brief analysis of situations where 
the velocity is constant and the system's behavior corre- 
sponds to that of a moving front, as well as others where 
the velocity of the bubble changes continuously. For the 
former case, and as mentioned in we recover a criti- 
cal velocity beyond which the front will no longer survive, 
see bold curve of As 5 ^ oo at the boundaries of the 
bubble, the curves will abruptly go to zero, rather than 
asymptotically as above. The analysis of this situation 
can be done by making a change of variables on Eq. Q . 
If we set ^ = x~ct we obtain the following equation which 
is valid in the moving frame accompanying the bubble. 



de 



(2) 



By measuring the maximum population density, 
of the front solution or the steady solution in a moving 
frame, we obtain a curve that will be useful in interpret- 
ing the following results. As the bubble accelerates in 
the moving frame, the system is not allowed to relax into 
a steady state, see the accelerating curve of FigQ] The 
time involved in the changes induced by the acceleration 
of the bubble compete with the relaxation time of the 
system. We will leave these results for now and recall 
them later to understand the following analysis. 



OSCILLATING BUBBLE 

Here, the habitable bubble of fixed width oscillates as 
a whole. We study two cases: sinusoidal velocity, and 
a velocity which is constant but suddenly changes direc- 
tion at a given distance away from the original position. 
All points within the amplitude of oscillation will thus 



be covered periodically for varying amounts of time. For 
the purpose of comparison, we will consider the velocity 
to be an average velocity for the sinusoidal case. The 
parameters involved are the average velocity of the bub- 
ble, V, the amplitude of oscillation characterized by the 
distance between the center of the bubble at the turning 
points. A, the length of the bubble, W, and the Fisher 
parameters: the diffusion coefficient, D, the growth rate, 
a, and the competition term, b{x,t). For the following 
analysis, b{x,t) is a step function that defines the shape 
of the bubble. We have examined the following aspects 
ofEq. CJ: 

• The role of each of the three parameters, v, A, and 
W, in the dynamics of the population. 

• How the two cases (sinusoidal velocity and constant 
velocity) differ from each other. 

• Qualitative characteristics and comparisons be- 
tween solutions for various parameter sets. 

• The evolution of the population as diffusion goes to 
zero. 

• The appearance of a curious dipping phenomenon. 



Numerical Results for D>0 

It is helpful to divide the behavior into two regimes 
which correspond to no-extinction and extinction. The 
parameters A, W, and the critical bubble width, Wc — 
TTy/{D/a) are sufficient to define these regimes. In 

the steady state, bubble widths smaller than the critical 
value are unable to support a population density. This 
feature now manifests itself in the amount of overlap be- 
tween the two extremes of an oscillating bubble. For 
A < W there will be an area that is permanently cov- 
ered by the bubble. When this covered area, from now 
on referred to as the overlap, is greater than the critical 
width, the population will always survive to some degree. 
This regime is bounded by < ^ < {W — Wc). 

In the other regime, extinction may occur when there is 
an overlap that is less than the critical value, (W—Wc) < 
A. The distance between the inner edges of the bubble 
when A > W will be called the separation, {A — W). 
Naturally, when there is a separation, no single location 
will receive continuous coverage by the bubble. 

For either regime, if the velocity is below the Fisher 
velocity, Vf, the population will be able to follow the 
oscillating bubble. As the velocity increases it becomes 
more difficult for the population density to keep up with 
the moving bubble (via diffusion and new growth) . Some 
of the population falls prey to the harsh environment 
which quickly reduces the density. If b{x, t) is infinite for 
|x|>W/2, the exposed population will go extinct, never 
to return. For realistic as well as numerical purposes. 
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FIG. 2: Temporal evolution of Umax where in (a) A = 8Wc, 
(b) A = 10 Wc, and W = 10 Wc in both. Each curve is char- 
acterized by one of seven sinusoidal bubble velocities, .laWc, 
laWc, SaWc, 5aWc, 7aWc, and lOaWc where Vf = aWc when 
D = aWc , as in the plots above. Higher velocities induce 
lower maximum population density values and thus drive a 
population closer to extinction as the amplitude of oscillation 
increases. 

we will consider b outside the bubble to be a large but 
finite number. Once outside of the bubble, the popula- 
tion decreases immensely yet it now retains the ability 
to regenerate. Due to the extended amount of time that 
either end of the bubble's trajectory is sheltered, peaks of 
increased population density form at these extremes. As 
the velocity is increased, these peaks become more pro- 
nounced with respect to the maximum population den- 
sity in the center of the bubble's trajectory. They also 
come more frequently since the bubble traverses the same 
period of oscillation more often, see in Fig. 12 and Fig. 
El The entire density profile oscillates spatio-temporally 
as the velocity increases and varies with the amount of 
overlap or separation. 

For the zero extinction regime, high velocities effec- 



tively wash out the oscillatory behavior and the peaks 
merge so that the length of permanent overlap then de- 
termines the population density. The behavior of the 
population in the extinction regime, {W — Wc) < A, is 
the same except for one crucial characteristic. Again, the 
peaks at the endpoints come closer together for higher 
velocities, but they never merge. The population dies 
out before a constant density forms in the center. As A 
increases, extinction occurs for smaller and smaller ve- 
locities. 

In the discussion above an average velocity of the bub- 
ble was assumed. We now make a distinction between 
two types of velocities and their relative affects on the 
population density. A bubble moving with sinusoidal ve- 
locity lingers for longer time over the end points of its 
trajectory and for shorter time in the center as compared 
to the constant velocity bubble. For a single location, 
a larger population density is permitted to grow if more 
time is spent under the bubble. The basic qualitative dif- 
ferences of the population density thus follow from these 
two aspects. Though the two velocities produce the same 
essential behavior, we find that there indeed exists a dif- 
ference between the two after long times. This will be 
addressed in the final section. 

To further understand the structural changes of the 
population, we have examined the maximum population 
density as it oscillates parametrically in time. The shape 
that is outlined can be roughly characterized by a bowtie, 
Fig. The densities of the maxima, minima, center, and 
the spatial extremes are the four features that comprise 
the qualitative differences between bowtie morphologies. 
As the average velocity increases these four points of 
comparison within the bowtie undergo changes in pro- 
portion. The maxima of the bowtie move closer to the 
center (recall that in the no-extinction regime these max- 
ima eventually merge). The minima move farther from 
the center. The spatial spread of the bowtie is smaller 
and the overall values of the population density decrease 
for all positions. The difference between the maxima and 
minima of the bowtie, Ai, (the amplitude of oscillation of 
the maximum population density) initially increases and 
then decreases as velocity increases. Fig. 0] depicts the 
abrupt appearance and then gradual decay of the bowtie 
as a function of velocity. The bowtie appears in both 
nondiffusive systems and diffusive systems. Bowties of 
the former develop amplitudes larger than those of the 
latter. Also, bowties are formed at lower velocities for 
the nondiffusive than for the diffusive system. 



Analytic Results for D=0 

We begin by considering the Fisher equation where 
the oscillation of the bubble is explicitly built into the 
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FIG. 3: A variety of bowties exhibit the affects for different 
values of v when a bubble oi A = 8Wc, W = lOWc and 
moves with sinusoidal velocity. Velocities from 
top to bottom range from aWc to lOaWc and increment by 
.5aWc where Vf = aWc when D = aWc- The trajectories 



are parametric in time such that u„ 
endpoints. 



is increasing at the 



FIG. 4: The amplitude of the bowtie is plotted against the 
average velocity (which here varies sinusoidally). The inset 
shows the average Umax after long time. The combination of 
the both show the validity of the analytic solution in the high 
velocity regime when the curves for D = aW^ and D = 
converge. 



argument of the competition parameter b{x,t) of Eq.lQ. 
--a u{z, tyb{z+ f v{s) ds) u{z, tf+D 



dt 



The transformation (j) = 1/u yields an equation of the 
form 



dt 



+ a(j){z,t) 



b(^z + J v{s) ds^ + D 



d^(t>{z,t) 
dz^ 



(j) \dz 



(3) 



By allowing Z? = 0, we recover the equation for a highly 
damped linear oscillator that is driven by an external 
force proportional to b{z + v{s)ds). The population 
will thus follow the periodic forcing function with a lag. 
This oscillator interpretation clarifies the structure of the 
bowtie. Each location grows periodically, with a tempo- 
ral lag, at the behest of the driving term. The locations 
at the edges of the bubble's spatial trajectory will have 
higher population densities because the driving oscilla- 
tor lingers for longer time at a lower value (less harsh 
conditions). 

An additional transformation 



v{s) ds, 



and letting D = leads to, 



The solution, as derived in |JJj, was found to be. 



(j)ix,t) = e 4>o{x- / v{s)ds) 
Jo 

+ j^e-^^'-^'nlx- v{s)ds\ dt' (5) 



and finally, u(x,t) explicitly as 



u(x, t) 



(6) 



uo{x - /p v{s) ds) 



dt 



dx 



e-^^ +uo{x- y^v{s)ds) /g e-° (*-*') 6(x-/J v{s) ds) dt' 

As the bubble's velocity greatly surpasses the Fisher ve- 
locity, the effect of diffusion becomes negligible. The dif- 
fusion length per time competes with the velocity of the 
bubble. Diffusion is thus limited by the time a single lo- 
cation is covered by the bubble and has little effect when 
velocity is high. A solution thus exists in the high veloc- 
ity limit, refer to Fig. 01 

When the velocity is relatively low, v < Vf, popu- 
lation migration occurs via diffusion and new growth. 
By letting Z) — > (correspondingly, Vf ^ and criti- 
cal overlaps 0) one might conclude that a population 
cannot maintain itself unless there is a bit of overlap. 
This is true when the velocity is such that growth can- 
not occur within the short period of coverage by the bub- 
ble. However, when these two rates (bubble velocity and 
aW) are comparable, oscillations in the population den- 
sity again arise in the diffusionless case. The bowties 
in these circumstances vaguely resemble those of before. 
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FIG. 5: Amplitude of the front oscillations for several veloci- 
ties of the bubble where D = aW^- 

When A < W the maximum population density is al- 
ways saturated; there are no means by which to leave a 
sheltered region without diffusion. When A > W, the 
bowtie limits to a J7 shape. 

BREATHING BUBBLE 

The breathing bubble corresponds to a situation when 
the width of a anchored bubble varies in time. As before, 
we consider several cases and for each of them a family 
of values for the involved parameters. These parameters 
are the velocity of the variation of the width, the mean 
width of the bubble, the amplitude of the breath, and 
again, the Fisher parameters of Eq. 1^. We are inter- 
ested in analyzing the relationship between the existence 
of critical values and the response of the population den- 
sity to a changing environment. We first take the case 
when the bubble width oscillates at a constant speed but 
always preserves a size above the critical one. The inter- 
esting aspect of this situation occurs at the boundaries 
of the bubble where population fronts are formed. The 
fronts can follow the movement of the bubble when the 
velocity of the breathing motion is under a critical value. 
Above this critical value, the population is confined to 
the minimum size of the bubble, which indicates that 
the population front can no longer follow the changes of 
the bubble. To show this we display in Fig. Elthe rela- 
tive amplitude of the oscillation of the population front, 
Af/A vs. the speed of the bubble breath, where A is 
the amplitude of breathing. We observe that at a critical 
value this amplitude goes to zero. That is, in spite of 
the changing size of the bubble, the population density 
remains bounded statically. 

Instead of a constant value for the velocity of breath- 
ing, we can take a fixed frequency for a sinusoidal behav- 



FIG. 6: Amplitude of the front oscillations for several values 
of the inverse of the frequency of breathing of the bubble 
where D = aW^. 

ior. The frequency is associated with the breathing mo- 
tion of the bubble boundaries. The results are analogous 
to those displayed in Fig. El No qualitative differences 
were found. 

When the minimum size attained by the bubble during 
the breathing is under the critical size, the situation is 
different. It is not enough to plot, as in previous exam- 
ples, the amplitude of the oscillation of the front because 
some variations in the maximum density of the popula- 
tion are also observed. To visually understand what is 
happening we again plot the temporal behaviour of max- 
imum value of the population, Umax as a function of the 
momentary width of the breathing bubble. This is what 
is shown in Fig. [3 where the trajectories are traversed 
clockwise. As observed in the plots, if the velocity is low 
enough, the size of the bubble will be under the critical 
size long enough to provoke the extinction of the popu- 
lation. On the contrary, for higher values of the velocity 
there is a range of values at which the population can sur- 
vive. The response of the population to changes is much 
slower than the dynamics of the environment. Each of 
the results obtained in this case are essentially the same 
as the corresponding results from the oscillating bubble. 
In particular, the results displayed in Figs. |3 and El can 
be associated with the discussion included in the previ- 
ous section, those included in Fig. Elare analogous to the 
bowtie effect. 

FURTHER FEATURES 

We now return to a surprising feature that arises only 
in the presence of constant bubble velocity and diffu- 
sion. Probing longer times reveals an overlying oscillation 
in the population density profile. For diffusive systems. 
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FIG. 7: Maximum value adopted by the population as a func- 
tion of the instantaneous width of the bubble. The veloc- 
ity of breathing is: a)0.1aW^c h)laWc c)10aWc d)40aWc, and 
D = aW!. 
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FIG. 8: After long times a periodic dip in all values of Umax 
is observed for constant velocities. The sampling above, (A, 
B, C, D, and E) corresponds to velocities {3.5aW, 4.5aW, 
5.5aW, 6.5aW, and 8.5aW, respectively) where «/ = aW 
because D = aW^. 

D > 0, and constant velocities, the entire bowtie dips 
down and then back up again in a periodic manner with 
constant ampUtude. For different velocities the bowtie 
dips and rises with different periods and amplitudes, see 
Fig. IHl Also, notice that as velocity increases, the fre- 
quency of the alternative dip doesn't simply increase as 
well. Nor does the amplitude of the dip appear to change 
with the velocity. Sinusoidal velocities do not produce 
this affect. Further analysis is necessary to completely 
understand the origin of this behavior and its relation to 
the two required conditions, diffusion and constant bub- 
ble velocity. 



CONCLUSIONS 

In the present work, we have numerically studied the 
behavior of a population whose dynamics are described 
by a modified Fisher equation with additional environ- 
mental variations. Analytic solutions were also consid- 
ered for simplified examples. The results presented here 
show that besides the expected behavior some surprising 
features arise. There exist two main parameters that 
characterize the evolution of a given population: the 
Fisher length and Fisher velocity. A given group of in- 
dividuals will not survive if the habitat size is smaller 
than the Fisher length. At the same time, the popula- 
tion will not be able to maintain saturation when the 
habitat moves with a velocity higher than the Fisher ve- 
locity. With these two considerations we have examined 
two types of periodically varying environments. Though 
extinction of the population was verified in conditions 
corresponding to the two situations mentioned above, it 
was also found that some nontrivial scenarios can also 
arise. Situations in which extinction was expected ex- 
hibited what we called the bowtie effect. An interplay 
between the time associated to the transient towards ex- 
tinction and the periods of the changing environments 
is the cause of unexpected behavior of the population. 
We have also verified that there is a correspondence be- 
tween the results obtained when considering oscillatory 
or breathing bubbles, though the situations are different. 
A diffusionless approximation helped to understand the 
origin of the observed phenomena in terms of oscillating 
perturbations. Perhaps the most unexpected result is 
that depicted in Fig|H| An explanation for these features 
is still in progress. 
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